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^ ■ We present results from a Monte Carlo simulation of non-compact lattice QED in 

O I 3 dimensions in which an explicit anisotropy k between x and y hopping terms has 

I been introduced into the action. Using a parameter set corresponding to broken chi- 

ji^! ral symmetry in the isotropic limit k = 1, we study the chiral condensate on 16^, 

CD . 20^ and 24^ lattices as K is varied, and fit the data to an equation of state which in- 

corporates anisotropic volume corrections. The value at which chiral symmetry is 
^ - apparently restored is strongly volume-dependent, suggesting that the transition may 

' be a crossover rather than a true phase transition. In addition we present results on 

16 lattices for the scalar meson propagator, and for the Landau gauge-fixed fermion 
propagator. The scalar mass approaches the pion mass at large k, consistent with 
chiral symmetry restoration, but the fermion remains massive at all values of k stud- 
ied, suggesting that strong infra-red fluctuations persist into the chirally symmetric 
regime. Implications for models of high-Tc superconductivity based on anisotropic 
QED3 are discussed. 

PACS: ll.lO.Kk, ll.15.Ha, 71.27.+a, 74.25.Dw 
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1 Introduction 



QED3, i.e. quantum electrodynamics restricted to two dimensions of space and one 
of time, has recently been the focus of some attention in the condensed matter com- 
munity, as various versions of it are examined as candidate effective models of high- 
temperature superconductivity in cuprate compounds. 

The particular model in which we are interested is that presented in |Tl |2] , which 
is supposed to model the passage between the antiferromagnetic SDW (spin-density 
wave) and superconducting dSC (the d denotes the superconducting order parame- 
ter has d-wave symmetry) phases at low temperature T as the doping fraction x is 
increased. QED3 is proposed, as reviewed below in Sec. 12.11 as an effective theory of 
the low-energy quasiparticle excitations in the neighbourhood of the 4 nodes in the 
gap function A(/c). Since the dispersion relation is linear at the nodes, the excitations 
can be reinterpreted as various components of a relativistic spinor field ^ with 4 spin 
and Nf = 2 flavour degrees of freedom. Interaction via a minimally-coupled abelian 
vector gauge potential field arises as a result of phase fluctuations of A; it can 
then be argued that is most naturally governed by an action of a Maxwell type 
Pl|2j, resulting in massless photon degrees of freedom which have an alternative inter- 
pretation as the Goldstone bosons associated with the condensation of dual vortices 

[3 HI- 

QED3 is a quantum field theory whose study has a long history (see j3] for a brief 
review). The main issue is chiral symmetry breaking (xSB), i.e. whether chiral sym- 
metry, the invariance of the action under independent global rotations of left- and 
right-handed helicity spinors, is spontaneously broken, signalled by a chiral conden- 
sate (^\E') 7^ 0. xSB implies dynamical mass generation, i.e. the physical fermion 
mass M may be much greater than the "bare" or Lagrangian mass m. This is be- 
lieved to depend sensitively on the number of fermion species Nf in the model; xSB 
is supposed to occur only for Nf less than some critical Nfc, whose precise value 
remains a goal of non-perturbative quantum field theory. 

In the condensed-matter context, the x^B order parameter can be mapped directly 
into the SDW one. If xSB does not occur (i.e. Nf > Nfc), then the resulting theory of 
light fermion degrees of freedom interacting with massless gauge degrees of freedom 
is proposed as a theory of the so-called "pseudogap" region of the cuprate phase 
diagram, characterised by spectral depletion in the immediate vicinity of the Fermi 
energy even in the absence of a well-defined quasiparticle peak. The main prediction 
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of the QED3 approach is thus that if Nf^ > 2 the dSC and SDW phases are connected 
in the T ^ hmit |2I , whereas if Nfc < 2 they are separated by a region of pseudogap 
phase pp. 

An important assumption in the above chain of reasoning is that results from con- 
tinuum QED3 in the isotropic hmit (as usually studied in quantum field theory) can be 
applied directly to the condensed-matter system, whose Lagrangian density (jH)) below 
has kinetic terms describing a single flavour with differing strengths or "velocities" 
in X- and ^/-directions as an artifact of the transformation to the relativistic spinor 
basis. This is in principle not a negligible effect; the velocity ratio or anisotropy k in 
real cuprates varies with x and can be as large as 7 at the onset of the dSC phase 
[Zj. Evidence in favour of applying predictions of the isotropic theory comes from a 
renormalisation-group analysis, which studied small anisotropy perturbations to the 
isotropic system and concluded that weak anisotropy is an irrelevant perturbation 
[Hlini- This result is then used to argue that the critical Nfc is a universal constant, 
independent of k, and hence that the various estimates of Nfc in the literature can 
be applied to the cuprate problem. 

It should be noted here that similar ideas regarding relativistic fermions (of- 
ten four-Fermi theories) have been also been discussed in the literature relating to 
graphene and similar compounds, both theoretically (for example, see [10111111121113) 
and experimentally (for example ^1]). However, this form has only an asymmetry 
between temporal and spatial directions, whereas in what follows we treat a more 
generally anisotropic system where all three Euclidean axes are distinguished. 

In our previous paper we made the first study of anisotropic QED3 using the 
methods of numerical lattice gauge theory, a non-perturbative technique with very 
different systematic approximations to continuum-based approaches. Using a large 
value of the fermion-photon coupling strength, we studied the xSB order parameter 
(xx) oil relatively modest 16^ spacetime lattice as the anisotropy k is increased 
from 1, and provided preliminary results that suggest the existence of a chiral sym- 
metry restoring phase transition at a critical value Kc, moreover, we found that the 
"renormalised" k/j - obtained by considering the spatial decay of correlations of pseu- 
doscalar meson or "pion" fermion - anti-fermion (//) bound states - obeyed > k, 
suggesting in contrast to that k is in fact a relevant parameter. Both observations 
suggest caution should be used in applying isotropic QED3 directly to cuprates. 

Several questions raised by are addressed in the current paper. Firstly, we 
wish to understand the nature of the chiral symmetry restoring transition, using the 
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traditional method of studying the transition on larger systems and applying a finite- 
volume scaling analysis. New results on 20^ and 24^ lattices, together with analyses 
assuming both isotropic and weakly anisotropic finite volume scaling are presented in 
Sec. El As we shall see, the data is best fitted by an ansatz which takes anisotropy 
into account, which suggests that the value of Kc in the thermodynamic limit may be 
considerably larger than the estimate — 4.5 of ^Hl- Next, in Sec. E] we have studied 
spectroscopy in another // channel with scalar, rather than pseudoscalar, quantum 
numbers. This is important for two reasons. Firstly, as the parity partner of the 
pseudoscalar the scalar should become degenerate with the pion at large k, giving 
further evidence for the restoration of chiral symmetry. Secondly, since the pion is 
the Goldstone boson associated with xSB in the 1ow-k phase, it is in some sense a 
"distinguished particle", and motivates us to check the renormalised anisotropy 
using a different channel. 

Finally, for the first time we present results for the fermion propagator {x{^)'X.{y)) ■ 
Since this is not a gauge invariant object, in order to obtain a non-zero result this has 
necessitated the implementation of a gauge fixing procedure, described in some detail 
in Sec. 15.11 Our motivation comes from the arguments of Tesanovic et al. ^ ITHj . 
suggesting that the massless quasiparticles of the pseudogap phase acquire a small, 
gauge dependent anomalous dimension due to their interaction with the statistical 
gauge field, which may explain non-standard scaling of transport coefficients such as 
resistivity and thermal conductivity in the pseudogap phase. From a numerical point 
of view this has proved easily the most demanding part of the project, requiring 
much computational effort to extract any kind of signal from the statistical noise 
inherent in the Monte Carlo method. Somewhat unexpectedly, we find evidence for 
the persistence of a djTiamically generated fermion mass in the high-K phase, despite 
the apparent restoration of chiral symmetry. A physical scenario consistent with these 
observations is discussed further in Sec. IHl 

2 Review of the Lattice Model 

2.1 QED3 as an effective theory of the pseudogap 

The mapping of the pseudogap region of the cuprate phase diagram onto QED3 
is derived in detail in |2j, and reviewed in language more accessible to particle 
physicists in fH]- Here we briefiy summarise, starting with the following Euclidean 
(imaginary time) action, also known as the Bogoliubov - deGennes model, for d-wave 
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quasiparticles in the dSC phase. 



A{k)cl{k, UJn)cl.^i-k, -UJn) - A\k)c„{k, Un)c^a{-k, -UJn)^ 
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a 
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where c\c are creation and annihilation operators for electrons with spin cr = ±1, 
UJn = (2n — l)7rT are the allowed Matsubara frequencies, the function is the energy 
of a free quasiparticle (which thus vanishes for k on the Fermi surface), and A(/c) is 
the gap function, which can be thought of as a self-consistent pairing field. Due to its 
(i-wave symmetry, A actually vanishes at two pairs of node momenta k = ±Ki, ±K2, 
with K1.K2 = 0. 

Linearising the latter functions around the nodes and defining the 4-spinor at 
the node pair i as 

u) = (c+(fc, 00), cl{-k, -00), c+{k - 00), cl{-k + -u)^ , (2) 

we may write the following effective action describing the behaviour of the system at 
low T 0: 

(3) 

where f3 = 1/T, k = vf/va (where vp and va are the Fermi and Gap velocities derived 
from the linearisation of and A(/c) respectively about the nodes) is the anisotropy, 
6 = y/vpVA, and the 4x4 traceless hermitian matrices 7^ obey {7^, -y^} = 25^^. The 
action (jHl) describes Nf = 2 flavours of relativistic fermion \1/ (sometimes known as 
'nodal fermions' in this context) interacting with an abelian gauge potential A^, which 
we will often refer to as the 'photon', via the covariant derivative = + iA^. 
The photon-fermion interaction models the effect of the phase fluctuations of the 
pairing field A: photon dynamics are governed by = {d[^Ay-^Y, and the coupling 
g (the analogue of 'electron charge' in textbook QED) is related to the diamagnetic 
susceptibility x via g ^ 111- 

The two velocities depend on the shape of the Fermi surface, and hence on the 
doping of the superconductor jHl El , implying that the same is true of k; at the onset 
of superconductivity at low T k may be as much as 0(7). 



5 



2.2 Lattice Model of Anisotropic QED3 

The formulation of isotropic QED3 on a spacetime lattice is described in detail in 
[T7] : in what follows we summarise the treatment of our anisotropic model given in 
[T3] . For flavours of staggered lattice fermion, the following is a QED3 action with 
an explicit spatial anisotropy: 

1 = 1 x,x' X,fJ,<U 

We define the fermion matrix M^^^/ as follows: 

1 ^ 

M^. = ^ ^m(^) [^x',x+fLU^fM - Sx',x-iiUl^] + m6f,^ (5) 

where is 

^t.{x) = Xf^Vi^ix) (6) 

and rif^{x) = where Xi = x, X2 = y and X3 = r, is the Kawomoto-Smit 

phase of the staggered fermion field. The physical lattice spacing is denoted by a. 
The are anisotropy factors, which we define like so: A^; = k~^, Xy = , Xt = 1. 
The 7] factors ensure that the action describes relativistic covariant fermions in the 
isotropic limit k = 1. 

Taking the photon- like degree of freedom 6f^{x) to exist on the link connecting site 
X to site X + fi, makes U^{x) = exp{ia9^{x)) in (jS]) the parallel transporter defining 
the gauge interaction with the fermions; we may define a non-compact gauge action 
via 

0,.(x) = l[A+^^.(x)-A+0,(x)]. (7) 

The dimensionless parameter /3 is given in terms of the QED coupling constant via 
(3 = 1/g'^a. It is convenient to work wherever possible in 'lattice units' such that 
a = 1. 

It is important that a distinction is made between our particular use of anisotropy, 
which treats it as a physical property of the system that can be observed and renor- 
malised through quantum corrections, and the more general use of anisotropic cutoffs 
in lattice field theory, wherein the anisotropics are controlled such that they disap- 
pear in the continuum limit, maintaining the Lorentz covariance of the theory. In 
this latter case, anisotropy is not physically observable. 
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2.3 The Simulation 



In Reference jlSj we simulated the dynamics of the lattice action ()4p5|l using a hybrid 
Monte Carlo algorithm on a 16^ lattice for k ranging from 1 to 10 and the bare 
mass m = 0.05, 0.01. The gauge coupling constant /? was held at a constant 
value 0.2 throughout - at this relatively strong coupling the system is in a state of 
spontaneously broken chiral symmetry at k = 1. The main results of jT^ are that the 
chiral condensate decreases with increasing n, consistent with a second order chiral 
symmetry restoring transition at = 4.35(2), and that the renormalised anisotropy 
kr obtained by comparison of pion correlators in x- and ?/-directions obeys 



implying that k is a relevant parameter. 

In the calculations presented in this paper, unless otherwise noted, the gauge 
configuration ensemble {6} used was generated using the same Hybrid Monte Carlo 
algorithm, running for around 1000 trajectories of mean length 1.0 on lattices with 
L = 16, and the gauge coupling set to the same value /? = 0.2. Even-odd partitioning 
was used; this allowed us to set = 1, giving us Nf = 2 in the continuum limit. 
Typical acceptances were 60 — 70% for m = 0.01, and 70 — 80% for other bare mass 
values. 

A novelty of this paper is that we have extended our study to a range of volumes: 
datasets for the 20'^ and 24^ lattices typically contain 700 and 600-700 trajectories 
per point with acceptances of 82 — 94% and 75 — 82% respectively. In the the studies 
of fermion propagation presented in Sec. El gauge-fixed configurations were generated 
on a 16^ lattice and consisted of ~30,000 trajectories per point with an acceptance 



We begin the presentation of our results with measurements of longitudinal suscep- 
tibility and the chiral condensate as L is varied. This is necessary to pin down 
the nature of the chiral symmetry restoring transition with more precision. Apart 
from the intrinsic theoretical interest, there are important phenomenological issues 
at stake. Firstly, it is important to know the value of the critical anisotropy Kc at 
which the transition takes place in the continuum and thermodynamic limits, since 
in principle this is a physically observable parameter in real cuprates [Oj. Secondly, 





rate of 79 - 87%. 



3 Susceptibilities and Finite Size Scaling 
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the order of the phase transition is important; were it either first-order or a crossover, 
then an immeasurably small but non-vanishing condensate may persist in the high-K 
"chirally restored phase" , meaning that antiferromagnetic order can survive the tran- 
sition As we shall see below, the results we have been able to obtain with our 
resources have not settled the issue unequivocally; it seems likely that a model of 
finite volume scaling which takes account of the anisotropy is required. 

3.1 Finite size scaling of the condensate 

Here we present the results of a preliminary study of the finite size scaling of the chiral 
condensate and longitudinal susceptibility at fermion mass m = 0.01, the smallest of 
the bare masses examined in ^21- We define the chiral condensate in terms of the 
trace of the inverse of the fermion matrix M: 



Note that eqn. (fTIH) includes diagrams which are both connected and disconnected in 
terms of fermion lines; both contributions were calculated. 

In the vicinity of the phase transition xi should peak at an anisotropy Kpeak which 
should tend towards the critical value in the thermodynamic limit. Examining 
the plot of the longitudinal susceptibility as the size of the lattice is varied (Figure 
HJ, we observe that the peak shifts to the right by an amount that decreases as the 
lattice size increases; this suggests that a second order transition might occur at a 
finite value of in the thermodynamic limit. Unexpectedly, however, the magnitude 
of the peak appears suppressed as the lattice size increases. This may have several 
possible causes: 

• The magnitude of the peak does increase, but the width of the peak as the 
lattice volume increases narrows such that it falls between the available data 
points and is not detected. The rounded shape of the curves suggests that this 
is unlikely. 

• This is not a second order phase transition; perhaps we're observing a crossover 
instead. If there is a crossover between the two phases, and not a genuine second- 
order phase transition, then a small chiral condensate is expected to persist in 



(9) 



and the longitudinal susceptibility in terms of its derivative. 




= -[((trM-i)(trM-i)) - {iiM~^Y - (tr(M-iM-i))]. 



(10) 
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9 mass(bare)=0.01, 16 cubed 
□ mass(bare)=0.01, 20 cubed 
^ mass(bare)=0.01, 24 cubed 

" 2 4 6 8 10 
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Figure 1: xi for various lattice sizes and m = 0.01. 

the high-K phase. Some analytic approaches predict dimensionless condensates 
P'^ixx) as small as O(10~^) [TH]. Attempts to rule out this possibility regarding 
the chiral phase transition in studies of isotropic QED3 with various Nf such 
as [T71 13] have not proven to be successful, and it is also likely to be as difficult 
in this case. 

• Our system has an anisotropic coupling between the gauge and fermion fields. 
The effects of this could be difficult to account for in the standard finite size 
scaling developed for phase transitions in isotropic systems; we should turn 
our attention to the scaling of anisotropic systems instead. In the statistical 
mechanics literature, one observes two models of this scaling: 

Weak anisotropy: In these systems, there exist different correlation 
lengths in different directions; these correlation lengths can be rescaled 
such that the system is effectively isotropic in the scaling region (^^1 and 
references therein, notably [201 1^ 122] )• We examine this possibility in 
detail below. 

Strong anisotropy: In these systems in addition to correlation lengths, 
the critical exponent u is different in different directions ([22, 23 and ref- 
erences therein). The scaling behaviour of these systems is very sensitive 
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to the shape of the lattice, and is difficult to treat with data generated on 
cubic lattices; it is mentioned here as an issue worthy of further investiga- 
tion. 

Since our data is restricted to that generated on a square lattice, we will examine 
only weak anisotropic scaling as compared to isotropic scaling. 

3.1.1 Isotropic scaling 

Firstly, we shall discuss the scaling of the system if it is assumed intrinsically isotropic. 
As argued in [21], we can use the scaling behaviour of the system as we vary L in 
order to determine how the finite volume affects the equation of state. We do this by 
treating the inverse linear size of the lattice, as an irrelevant scaling field and 
use the following as our ansatz, where k = [k, — n^): 

m = B{xxY + A{k + CL-^'''){xxY- (H) 

Here (5, v and j3 = {5 — p)^^ have their usual meanings as critical indices describing 
a continuous phase transition. 

3.1.2 Weakly anisotropic scaling 

In this case, we wish to account for the distortion of the correlation lengths of the 
system along the x and y axes by k 7^ 1. Finite size effects enter into the scaling 
whenever 

» L,, (12) 

where is the correlation length in the direction /i and is the length of the lattice 
in that direction. We introduce three irrelevant scaling fields: Ll^^ and L^^ ^ 

defining them in terms of L, the number of lattice spacings along one dimension of 
the system, by rescaling [13 EOl EB E2! such that 

Ci = C2 = (13) 

Since (to a first approximation) = ^/K = ^/k^ and SJf = ^ (where ^ is the 

correlation length of the isotropic system), this gives 
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This rescaling of the correlation lengths is equivalent to resizing the lattice thus (from 
consideration of (|T^ ): 

Li = \/kL, L2 = a/k ^L, L3 = L. (15) 

So, we can write: 



+D^\ +E 



L J \ L J \ \fKL ^ 

R{k-C,D,E)L^I\ (16) 



This motivates the replacement 

C 



R{k;C,D,E)L-^/'' (17) 



in (|lip . which we may then use to study the scaling if weak anisotropy is assumed. 
3.1.3 Results and Discussion 

We should note that the above equations are only good descriptions of the behaviour 
of the system near to a continuous phase transition. We have attempted fits to the 
finite-volume equation of state ()11|) using data from 16^, 20^ and 24^ with m = 0.01. 
To ensure stability of the fit we found that it was also necessary to include the 
m = 0.02 data for the 16^ lattice, presented in ^H], giving 34 data points in all. 

In addition, in order to increase the tractability of our fits, we have made use of 
the following hyperscaling relation (with dimensionality set to 3): 

, (^ + 1) 



3(^ - p) 



which reduces the number of free parameters in our fit to six assuming isotropic 
scaling and eight assuming weakly anisotropic. 

Results from fitting the chiral condensate data to (fTTj) are shown in Table H The 
daggered quantities were obtained through the following relations: 



11 



Quantity 


Isotropic 


Weakly anisotropic 


A 


.0393(8) 


.0111(7) 


B 


1.28(8) 


1.02(6) 


C 


368(16) 


-755(338) 


D 


— 


527(78) 


E 




-716(447) 




( .DD(0 j 


iz.6(b ) 


5 


3.40(6) 


3.33(6) 


P 


.991(7) 


1.01(1) 




.41(1) 


.433(3) 




.363(6) 


.386(7) 




.61(2) 


.62(2) 


d.o.f. 


162 


6 



Table 1: Equation of state fit results, allowing for finite size scaling. Daggered values are 
calculated from hyperscaling relations (see main text). 




Figure 2: The corrected plot of the chiral condensate and the equation of state fits on a 
fixed volume 16^. 



The equation of state fits found on 16^ are plotted in Fig. |21 ^; for comparison 
the new equations of state, together with the fitted data and the extrapolation to 
the chiral limit m ^ 0, are plotted in Fig. 01 The following features are perhaps 

^in fact, the original figure shown in |15j had incorrect curves, and Hc was not located correctly, 
although the values of the critical exponents given were correct, and the conclusions of that paper 
remain unaffected. 
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the most intriguing. The critical indices are compatible for both fits - however, the 
value of Kc is not only different from the value Kc = 4.35(2) derived from fitting to the 
16^ data alone but is significantly different between the two forms of the finite 
scaling fit. This suggests not only that finite size effects play a significant role in the 
behaviour of this system, but that the effects of the anisotropy should be taken into 
consideration in future studies of the system. It is also worth noting that the value 
of the x^/d-o.f. is significantly better for the anisotropic scaling. Note also that for 
the weakly anisotropic fit, the sign of the coefficient of L2, D is different from those 
of Li and L3, C and E. This may reflect the expectation following (fT3j) that L 
over a much wider region of k than is the case for the other two directions. Contra 
the fitting results of [15^, which were confined to a single lattice size, p ~ 1.00 with 
both equations; however, it is plausible that this is due to an insufficient spread of 
mass values in the data set. 

Whether the actual value of the weakly anisotropic Kc is in fact 12.3(6) seems 
doubtful; we must note that the extrapolation is well outside the region of k, for 
which we have any data. An interesting possibility is that it could also indicate that 
there is no phase transition and that the fit could be attempting to compensate for 
its absence by giving it a value in the unexplored region. If this behaviour were to 
persist for a more extensive data set, this hypothesis could be validated. 
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4 Scalar Sector 



The scalar meson is the parity partner of the pseudoscalar pion bound state studied 
in In a phase with broken chiral symmetry the pion is a Goldstone boson, and 
hence is much hghter than the scalar. One signal for restoration of chiral symmetry 
is the recovery of degeneracy between scalar and pseudoscalar in the m — > limit. 
The propagator of the scalar is defined in terms of the fermion fields as follows: 

C^,{x,) = Y,Y.^XXmx{x)). (20) 

Due to the nature of the flavour structure of staggered lattice fermions, propa- 
gation in this channel is prone to mixing with low mass bound states with different 
spin quantum numbers Where this contamination is significant, the propagator 
takes on a sawtooth shape, and we must thus fit a four parameter function, such as 
that in ()22j) below, so that we can distinguish propagation in the channel of interest. 

In the following we distinguish between propagation in the Euclidean time direc- 
tion r, yielding information on the excitation spectrum in the channel in question, 
eg. the bound-state mass, and propagation in the spatial directions x, where the 
corresponding quantity is the inverse screening length. Of course, in an isotropic 
system the two cases are equivalent in the infinite volume limit. 

4.1 Temporal propagator 

Least squares fitting of the function 

C^,,{x^) = A(e-'"--"- + e-'"-''^^''-"^)), (21) 

(with n chosen to be r) to data from 16'^ lattices proved to be difficult within the 
chirally broken phase - the propagator data was exceedingly noisy, and care had to 
be taken in order to isolate the ground state signal from the excited states - but 
as the values moved into the chirally restored region the procedure became easier to 
perform. The results are listed in Table El and plotted in the graph El alongside the 
pion masses of for each bare mass at k = 10.00. 

It can be seen from the figure that there are two regimes of scalar behaviour; 
below Kc, where fitting is quite difficult, is more or less constant as k, increases 
(if we go by the m = 0.01 data and ignore the outlier at k = 3.00) up to k ~ 5 (ie. 

estimated on 16^), whereupon we find that nia begins to converge with rrin 
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Figure 4: Scalar masses in the r direction. Straight lines represent the pion masses at 
K = 10.00, taken from JH]- n^- = 0.05 values are omitted due to the size of their error bars. 



as K increases into the chirally restored region. Figure El shows this in more detail for 
m = 0.01. 

We should point out that the jump in the value of rricrT at k = 3.00 and m = 
0.01 (and likely that of rrifjy at k = 4.00, m = 0.05, see below) is likely to be due 
to the frequent occurrence of abnormally small eigenvalues of the Dirac operator 
that overlapped with the meson source during our measurement of the propagator, 
similar to that seen in the Thirring model simulations of [2^]. The overall trend is 
consistent with the scalar becoming degenerate with the pion at large k, consistent 
with manifest chiral symmetry. There is thus no evidence for persistence of chiral 
symmetry breaking at large n from the light meson spectrum. 



4.2 Spatial propagators 

The spatial scalar masses were also obtained by least squares fitting to the propagator. 
For the x-direction (Table El) and the k = 1.00 y-direction correlation functions, 
we used the fit function (jSH), and selected the fit window so as to exclude higher 
mass states. For k > 1.00, the y-direction correlation function exhibits a saw-tooth 
behaviour, motivating the following fit: 

C^y = A(e-™^«^ + e-"''^y^^y'y'>) + {-iyB{e-^'^y + e-*^(^^~^)), (22) 
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Figure 5: Scalar and pion masses (from [^) in the r direction for m = 0.01. 

which proved acceptable across the full range of k if a fixed fitting window of spaces- 
lices 1-15 was used. As with the temporal scalar propagators, those in the chirally 
symmetric phase were easier to fit than those in the chirally broken phase. It is also 
worth noting that in the symmetric phase the value of the correction mass M was 
often consistent with zero for m = 0.01 and m = 0.03. 

We have plotted m„ ^ against k in Figure El and mo- y against k in Figure 
The trends previously observed for pions in ^H] are repeated here: the value of max 
increases with n. In addition, as we have seen with maT, it appears that there 
is convergence between the rricrx and the rriTrx values, most notably for m = 0.01 
(compare m„r^ = .32(4) and m^rx = -211(1) at k = 1.00 with rriax = 2.62(2) and 
m.,rx = 2.58(2) at K = 10.00). Just as in Fig. 01 there is a change in behaviour around 
K ~ 5, suggesting a change in behaviour as the scalar masses begin to converge on 
the pion masses within the chirally symmetric phase, once again consistent with the 
pion and scalar being parity partners. 

In the case of the ^/-direction masses, there is no pion data within the chirally 
symmetric phase with which to compare our results. The quality of the scalar data 
is also not that good, for the reasons mentioned above. It is less clear whether the 
change in the behaviour between phases is present here; it is likely to be quite small, 
and in any case the errors easily obscure it. Our only conclusion is, then, that the 

16 
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m 


K 




xVd.o.f. 


fit window 


0.01 


1.00 


0.43(6) 


1.991 


2-14 




2.00 


0.40(8) 


1.882 


2-14 




3.00 


1.04(2) 


1.323 


1-15 




4.00 


0.47(3) 


1.045 


1-15 




5.00 


0.43(3) 


1.144 


5-11 




6.00 


0.55(1) 


1.119 


1-15 




7.00 


0.65(1) 


0.339 


1-15 




10.00 


0.91(2) 


0.984 


3-13 


0.03 


1.00 


0.8(2) 


1.177 


2-14 




2.00 


0.6(1) 


1.089 


2-14 




3.00 


0.7(2) 


0.923 


2-14 




4.00 


0.8(1) 


0.397 


2-14 




5.00 


0.75(4) 


1.351 


2-14 




6.00 


0.78(4) 


0.983 


3-13 




7.00 


0.82(2) 


0.505 


1-15 




10.00 


1.07(2) 


1.273 


1-15 


0.05 


1.00 


1.0(2) 


2.519 


2-14 




2.00 


1.2(7) 


1.036 


2-14 




3.00 


0.6(1) 


1.806 


2-14 




4.00 


1.0(2) 


0.745 


2-14 




5.00 


1.4(1) 


0.739 


2-14 




6.00 


0.95(5) 


1.376 


2-14 




7.00 


1.05(5) 


0.635 


2-14 




10.00 


1.16(2) 


0.662 


1-15 



Table 2: Scalar masses niscT in the r direction on a 16'^ lattice, for various masses. 

value of TJifj y decreases as we increase k; this is reinforced by the behaviour of the 
geometric mean of nifjy and ^ above k ^ 5 ~ (Fig. EI); while there is a slight 
increase as we approach Kc (though the large error bars make it difficult to determine 
to what extent this is a genuine effect), above it, the geometric mean appears to 
remain fairly constant, which implies that m^y is decreasing, as mg-x is increasing in 
this region. 

4.3 Renormalised Anisotropy 

Taking k/jo- to be the ratio of m„.j. and mg-y (Figure El), we find that it is relevant 
above /tc (more so, in fact, than for pions; Cf. Fig. 8 of ^Sl); however, in contrast to 
the pion case there appears to be a clear mass dependence below k^. It is difficult to 
tell whether this is a real effect or merely an artefact of the propagator fitting. On 
the assumption that the behaviour for m = 0.01 is more or less linear, we fitted the 
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m 


K 


rrisx 


9/1 £ 

X /d.o.f 


fit window 


0.01 


1.00 


0.32(4) 


1.052 


2-14 




2.00 


1.0(2) 


0.856 


1-15 




3.00 


2(1) 


0.874 


1-15 




4.00 


1.04(8) 


0.645 


1-15 




5.00 


1.28(3) 


1.215 


1-15 




6.00 


1.53(2) 


1.153 


1-15 




7.00 


1.83(2) 


0.734 


1-15 




10.00 


2.62(2) 


0.55 


1-15 


0.03 


1.00 


0.7(2) 


1.643 


2-14 




2.00 


1.2(2) 


1.495 


1-15 




3.00 


1.7(7) 


1.021 


1-15 




4.00 


1.8(3) 


0.577 


1-15 




5.00 


2.3(2) 


1.09 


1-15 




6.00 


2.2(1) 


1.596 


1-15 




7.00 


2.3(1) 


0.984 


2-14 




10.00 


2.80(3) 


0.853 


1-15 


0.05 


1.00 


0.7(1) 


2.485 


2-14 




2.00 


1(1) 


0.47 


2-14 




3.00 


2.2(4) 


0.981 


1-15 




4.00 


3(3) 


0.967 


1-15 




5.00 


2.5(4) 


0.842 


1-15 




6.00 


2.6(2) 


0.637 


1-15 




7.00 


2.8(2) 


0.906 


1-15 




10.00 


3.20(7) 


1.029 


1-15 



Table 3: Effective scalar mass x in the x direction. 



data for 1.00 < k < 7.00 to 



n, ^ (23) 



and acquired = 2.8(1), with x^/d-o.f. = 1.77. This appears slightly larger than 
Rt, ~ 2.1 [in], suggesting that the behaviour of scalar particles is affected by the 
anisotropy to a greater extent than that of the pions. 

More data is needed before we can make definitive statements. In addition, it 
should be noted that we can't rule out the existence of a change in behaviour around 
Kc for the values of krt^ - the pion data of ^Hl doesn't extend far enough. Based on 
the parity partnership of pions and scalars, we propose the following hypothesis: that 
in the chirally restored phase, the magnitude of k^jo- will gradually approach that of 
krt,. The required data would best be generated on considerably larger lattices, with 
better statistics and perhaps with improved operators in order to avoid some of the 
issues with the data examined here. 
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m 


K 


Ulfjy 


2/7 £ 

X /d.o.f 


fit window 


0.01 


1.00 


0.41(6) 


0.56 


2-14 




2.00 


0.25(5) 


1.469 


1-15 




3.00 


0.2(1) 


1.923 


1-15 




4.00 


0.13(2) 


1.548 


1-15 




5.00 


0.09(2) 


2.142 


1-15 




6.00 


0.1(2) 


1.227 


1-15 




7.00 


0.1(5) 


0.608 


1-15 




10.00 


00(11) 


1.735 


1-15 


0.03 


1.00 


0.8(2) 


0.721 


2-14 




2.00 


0.59(7) 


1.732 


1-15 




3.00 


0.7(3) 


1.076 


1-15 




4.00 


0.50(6) 


7.661 


1-15 




5.00 


0.15(4) 


1.015 


1-15 




6.00 


0.12(9) 


1.351 


1-15 




7.00 


0.10(8) 


0.798 


1-15 




10.00 


0.1(2) 


1.55 


1-15 


.05 


1.00 


0.9(2) 


0.676 


2-14 




2.00 


2.9(6) 


8.468 


1-15 




3.00 


0.8(3) 


4.308 


1-15 




4.00 


0.28(7) 


1.315 


1-15 




5.00 


0.19(5) 


2.097 


1-15 




6.00 


0.15(4) 


1.636 


1-15 




7.00 


0.12(5) 


1.818 


1-15 




10.00 


0.1(3) 


1.87 


1-15 



Table 4: Effective scalar mass rricry in the y direction. 



5 Fermion Sector 

In this section we report for the first time on studies of the fermion propagator 
{x{^)x{y)) ■ Large non-perturbative corrections to this Green function in the chiral 
limit m — > have been proposed as an explanation of non-Fermi liquid behaviour 
in the non-superconducting region of the cuprate phase diagram P]. An important 
challenge, both technical and conceptual, which must be faced is that the fermion 
propagator in QED is not a gauge invariant object, and can only by calculated, 
either analytically or numerically, if a gauge-fixing procedure is specified j^H]- The 
dependence of the results on the choice of gauge is a thorny issue j27| |2H1 CHI 123 ; here 
we will content ourselves with specifying Landau gauge, ie. d^A^j, = in continuum 
notation (implying that only transverse degrees of freedom are retained in the photon 
propagator), and performing a fully non-perturbative calculation on a 16'^ lattice. In 
what follows we will first devote some considerable attention on the technicalities of 
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Figure 6: The scalar screening mass in the x-direction, tUsx, on a 16'^ lattice. The straight 
lines represent rriTrx at a value of k = 10.00. 



fixing an unambiguous gauge for lattice gauge fields ?7^, and then report our results 
for the fermion propagator. Our strategy in this exploratory study is to calculate the 
physical (ie. renormalised) fermion mass ruf for fixed bare mass m as a function of 
the anisotropy parameter k. Apart from the fact that this is the simplest quantity 
to extract (by fitting to a decaying exponential), there is the theoretical motivation 
that rrif, given by the position of a pole in the complex fc-plane, is gauge-invariant, at 
least to all orders in perturbation theory. As previously, we will distinguish between 
propagation in temporal and spatial directions. 



5.1 Gauge Fixing 

In order for the measurement of a gauge variant quantity such as the fermion propa- 
gator to be performed, we must impose a gauge condition which selects a unique set 
of gauge configurations from the infinite number of copies generated by local gauge 
transformations of the form (in this subsection we will denote the lattice site by a 
suffix) 

0^.^e'^, = e^x + d^ax, (24) 
where on a lattice finite difference operators are defined: 

dfifx fx+ji fxi d^fx fx fx—jij (25) 
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Figure 7: The scalar screening mass in the y-direction, rrisy, on a 16^ lattice. We omit the 
anomalous m = 0.05, k = 2.00 value since its inclusion obscures the general trend of the 
data, and the m = 0.01, k = 10.00 value due to the size of its error bars. 

and ax is any scalar function defined on the lattice sites. 

For this study, we shall impose a latticised form of the Landau gauge condition 



(26) 



which is the extremum of 

^"[^] = EE(^P'' (27) 

X /J=l 

corresponding to the following functional in terms of continuum gauge fields: 

F[A] = J d^xA^{x)Af'{x). (28) 

In order to proceed, modifications will have to be made to this minimal gauge 
condition (henceforth referred to as mLandau gauge). This is because it suffers from 
the so-called Gribov ambiguity pUj . 

5.1.1 The Gribov problem in QED3 

When gauge fixing is performed non-perturbatively, it may always not be possible 
to guarantee that there is a unique minimum of the functional F[6]. In numerical 
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Figure 8: The geometric mean of the scalar screening masses, y/rus yfUgx-, on a 16'^ lattice. 



simulations, this can lead to a distortion of the results due to the underlying am- 
biguity jnij. The problem is normally associated with non-Abelian gauge fields in 
the continuum; however, it exists for Abelian fields on the lattice due to the toroidal 
boundary conditions which give rise to zero modes which cannot be removed by 
local gauge transformations and is especially acute for compact (cQEDs) formulations 
of the gauge fields as it allows for the existence of topological defects (such as double 
Dirac strings in 2 + 1 dimensions or double Dirac sheets in 3 + 1 dimensions) whose 
creation or annihilation leaves the action unchanged [5^ . 

Since we make use of a non-compact formulation of QED3 (ncQEDs ) in this study, 
it seems that the only problem we might have to deal with is the former. The modified 
iterative Landau gauge (miLandau gauge) jHH ESJ EH] has often been used in order to 
deal with the problems due to the existence of zero-modes created by the boundary 
conditions of the lattice; however, it has not (as far as we are aware) been checked 
that there are any other sources of Gribov copies in this gauge. So, in what remains of 
this section, we shall describe miLandau gauge and present results that demonstrate 
that it does deal with the problem effectively, at least for the values of the parameters 
simulated in this paper. 
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Figure 9: The renormalised k with respect to scalars, KRa, on a 16^ lattice, together with 
a linear fit for m = 0.01 (the dashed line is the quadratic fit, the filled line the linear). 
K = 10.00 has been omitted due to the size of the error bars. 



5.1.2 The miLandau gauge for ncQEDs 

Firstly, we note that on the lattice we cannot rotate 6"^ i-^ + a^, where is an 
arbitrary constant vector field, if we wish to preserve the gauge invariance of the 
Polyakov and Wilson lines (defined to be products of the parallel transporters U^^ 
along contours which are closed by periodic boundary conditions in the temporal and 
spatial directions, respectively). Instead, the form of the allowed gauge rotations is 
restricted to = where n is an arbitrary integer. 

Using 9fj_ = y J2x ^fj-x as the value of a constant background field (our zero- mode) 
we should expect the gauge degrees of freedom remaining after the mLandau gauge 
is fixed to vanish if we rotate 

0;^9; + ^ (29) 

such that <en< j^^. 

• We fix mLandau gauge using a steepest descent algorithm |HH] : 

similar prescription, the Zero-Momentum Landau gauge sets — 0. The difference 
between this and miLandau gauge in the thermodynamic limit (L^ oo) should be minimal j35| . 
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o Given a gauge configuration {9}, for eacli site we calculate tlie value of 

o If '^^Gx < R, where R the floating point value 10"^, we terminate the 
algorithm here. Otherwise, we continue. 

o We rotate 6*^3; i— 9f^{x) — d^Xx on every link of every lattice site, where 
Xx = V^x, and is a tunable parameter (here set to a value of 0.2), used 
to optimise convergence. 

o We repeat until the halting criterion is fulfilled. 

• Once mLandau gauge fixing is complete, we calculate for jj, = 1. 

o Add 1^ to each 9 ax until < 9^ < ^. 

o Subtract |^ to each 9u.x until -f- < < t-- 

• Otherwise, leave each 9^^ unchanged. 

• Repeat the above for the remaining directions /i = 2, 3. 

5.1.3 A test of this prescription in ncQEDs 

We wish to check that miLandau gauge removes Gribov copies from our measure- 
ments, by testing the effects of imposing miLandau gauge on randomly generated 
gauge copies of a set of gauge configurations [SHI HO] • The results were generated for 
K, = 1.00 and k, = 10.00 on a 16^ lattice for (3 = 0.2 and m = 0.03, the extreme values 
of the range at which we wish to measure the propagator. 200 mother configurations 
were generated, and for each mother we created three 500 configuration ensembles 
corresponding to one of the following random gauge transformations: 

Group A: For all x and fi, 9^^ ^ 9^^ — d^^ax where ax is a random number 
between —9 and 9. 

Group B: For all x and /x, 9^x ^ 9^x + ^/if^; where is a random integer - 
either 1, or —1. 

Group C: We perform both the transformation performed on Group A and that 
performed on Group B. 
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During the gauge fixing of eacli configuration, we monitored both the average 
value, Fav of the gauge fixing functional (fTTjl at each site and the value of the the 
function L = yJ2x^x, with Gx defined in ^5.1.21 for each iteration of the fixing. The 
behaviour of these parameters for a typical configuration is shown in Figure ITUl 
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Figure 10: Typical behaviour of Fav (a) and L (b) during gauge fixing. The plateau towards 
the middle of (a) corresponds to the area where (b) is converging on zero (that is to say, 
approaching mLandau gauge); the fall off beyond the hundred and fourth interaction is due 
to the imposition of full miLandau gauge. The first few points of both plots have been 
omitted so that this behaviour is visible. 

We may also define a 'variance' dF |37J, which measures the difference in the 
minimised values of the gauge functional -Fmm in a particular ensemble of a mother 
and associated daughter copies: 



dF = maxij[F„ 



mm ]] 



(30) 



with i,j = 1, . . . ,N where N is the number of daughter configurations in the ensemble. 
If there are no Gribov copies present, this quantity should be zero (more realistically, 
in a numerical simulation we expect it to be of the order of the residual, 10~^), 
otherwise we expect a large value. 

The results of our simulations at k = 1.00 and k = 10.00 with respect to Gribov 
copies were identical; we display figures for the former case, but our comments should 
be interpreted as generalising over both values of k. Figures [TT] and [T^l plot F^m and 
dF for the ensembles generated using each procedure. 

Group A: Here we find that while the value of -Fmm is appreciably different in 
miLandau gauge from that in mLandau gauge (indicating that zero 
modes exist and have been gauge rotated away in the former), dF is 
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Figure 11: Plots of F^m at k = 1.00 for Groups A, B and C. Fmin for B and C are identical 
to within 10~^. 




Figure 12: Plots of the variance dF at k = 1.00 for Groups A, B and C. dF for B and C 
are identical to within 10""^. 

of the order of 10~^, suggesting that the random gauge transformations 
used here do not usually generate Gribov copies. 

Group B: Unlike the above case, here we can see that there are in fact Gribov 
copies in the mLandau gauge: dF is between 4 and 5 for k, = 1.00 and 5 
and 6 for k = 10.00. However, this is not the case for miLandau gauge. 
Here, as before, dF in miLandau gauge is of the order of 10"'', and 
thus we can conclude that it rids us of the Gribov copies introduced 
by the random gauge transformation. 
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Group C: Here, the crucial observation is that (to within 10"'''), the results for 
these random gauge transformations are identical to those of Group 
B. Indeed, it would be worrying were it otherwise; the effects of the 
two sets of transformations should be additive, so one would expect 
only Group B's transformations to have any effect. 

It is clear from the above that only shifts in the constant background field appear 
to contribute to gauge copies, and these are readily dealt with through the addition of 
further constraints to the minimal gauge fixing condition, via the choice of miLandau 
gauge. This stands in strong contrast to the case of CQED3, where the compact 
Wilson gauge action allows for the existence of additional topological defects [HH! 
which are also solutions of the equations of motion and therefore are Gribov copies. 

This 'desert landscape' with respect to Gribov copies is not a disappointment - 
in fact, it is precisely the situation desired; one can be sure that the gauge has been 
fixed as in an unambiguous fashion. 

5.2 The Fermion Propagator 

In this section, we present measurements of the fermion propagator in the temporal 
and spatial directions: 

CM = 5^ $^(x(0)x(a:)), 

_ \ '^Vv Even numbered slice . ^ 

[ ly^ + /t Odd numbered slice 

where the sum on x only includes sites which are displaced from the origin by an even 
number of lattice spacings in each of the two transverse directions [231|2nj. We have 
also imposed noncompact miLandau gauge using the procedure outlined in §5.1.21 
As mentioned previously, the calculation of Cf^ required the generation of around 
30,000 trajectories of mean length 1.0 per K-point in order to extract a signal from 
the considerable noise; for a dynamical fermion simulation this amounts to a large 
effort, requiring between one and three weeks per point to complete. Because of this 
difficulty, the error bars of our measurements remain sizable. 
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5.2.1 Temporal propagator 

We extracted the fermion mass rrifr in the temporal direction from the propagator 
data via the function 

Cfrir) = Aie-'^f-^ - (-l)^e-'"/-(^--")) (32) 
using correlated least-squares fitting; the results are recorded in Table Eland Figure 

M 

First we should first examine Figure ^1 which shows examples of fermion propa- 
gators in the chirally restored phase k > 0.5. The following should be noted: 

• The central area of each propagator is fairly fiat, with large error bars (true 
of both phases). In this region the signal is overwhelmed by the noise and is 
consistent with zero. The size of the window containing data points exhibiting 
this behaviour decreased as the number of configurations in the sample was 
increased, suggesting that the cause is insufficient statistics. Because of this, it 
proved necessary to use fitting windows that are wider than the noisy region 
in order to extract a mass from the propagator. As in previous studies of ele- 
mentary fermion propagation Plj, no indication of contamination from excited 
states was seen within those windows. 

• Figure El also illustrates an interesting feature of the fermion propagators for 
K > 5: the onset of a sawtooth-type behaviour visible in the logarithmic plots 
that, although relatively small, grows more pronounced with increasing k. Since 
it is hard to distinguish it from noise, we performed fits of (j!?^ to i) all of 
the timeslices and ii) to only the odd numbered timeslices for the propagators 
exhibiting this behaviour. Ideally, a four-parameter fit is preferred to ii), but 
these proved to be unstable. 

The lines of best fit for both i) and ii) are included in Figure El for purposes of 
comparison, and the masses extracted are included in Table El and Figure IT^ 

It is worth discussing the origin of the sawtooth behaviour. The chiral symmetry 
preserved by the lattice model ()4f5j) in the limit m — >^ is the U(l) rotation 

X{x) ^ exp{i(3e{x))xix); x{x) ^ exp{i(3e{x))xix); (33) 

where the phase e{x) = (— i)^'i+^2+x3 (iigtinguishes between even (e) and odd (o) 
sites. In the chiral limit the only non-vanishing entries of the fermion propagator 
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Figure 13: Comparison of all timeslice and odd timeslice fits on (a) a linear scale, and (b) 
logarithmic scales for k of 6.00, 7.00 and 10.00 in descending order The errorbars represent 
unbinned, raw, statistical errors. 



matrix are M~^^ and M^^; for small but non-zero m it should still be the case in 
the chirally-symmetric phase that |M^-^|, iMg^-*^! ^ l^e'e^l; l^t^^l- ^^e timeslice 
correlator Cf^{x^) defined in (pTH) this implies that the signal should be much larger 
if is odd. Figure IT^ shows that the sawtooth behaviour of the curve is not especially 
pronounced, and it is at present unclear to what extent the phenomenon is connected 
with the restoration of chiral symmetry. 
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K 


m 


xVd.o.f. 


fit window 




1.00 


1.00(2) 


1.753 


2-14 




3.00 


1.17(2) 


1.696 


1-15 


fit 


4.00 


1.33(7) 


0.919 


2-14 


all timeslices 


5.00 






i-io 




6.00 


1.5(2) 


1.211 


2-14 




7.00 


1.3(2) 


1.031 


2-14 




10.00 


1.7(5) 


1.015 


2-14 


fit only 


6.00 


1.58(9) 


0.546 


1-15 


odd timeslices 


7.00 


1.6(1) 


1.228 


1-15 




10.00 


1.8(2) 


0.863 


1-15 



Table 5: Fermion masses m/7- in the r direction on a 16^ lattice, with m = 0.03. 



^ mass{bai"e)=0.03, fits all slices 

O mass{bai"e)=0.03, fits only odd slices 
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Figure 14: The renormalised fermion mass, mj, measured in Landau gauge on a 16'^ lattice 
for m = 0.03. 

Fig. 1141 shows that for n < 5.00, m^^ increases with k. The behaviour above k, = 5 
depends on the the type of fit - for fit ii) we see that the behaviour shows a non-zero 
mass in the region, which is more or less constant . Fit i) also shows the existence 
of a non-zero mass, but with more noise, possibly since it does not account for the 
sawtooth behaviour. 

Regardless of the method chosen for the fitting of the propagators, there is a clearly 
a non-zero dynamically generated fermion mass in the chirally restored phase. This 
is unexpected - dynamical mass generation usually implies (xx) 0^ chiral sym- 
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metry restoration usually implies massless fermions (Cf. Figs. 14 and 17 of Ref. |24j . 
illustrating light fermion propagation on a 16^ lattice in the chirally-symmetric phase 
of the 3d Thirring model). This seems to indicate that we are observing an unusual 
kind of chiral symmetry restoration; we shall return to this issue in due course. 

5.2.2 Spatial propagators 

We fit the spatial propagators to the following function: 

Cff,{x^) = A(e-'™^>^^ + (-l)^Me-™/M(iM-^M)) . (34) 

the change in sign compared to eqn. lp?^ being due to the use of periodic boundary 
conditions for the fermion fields in spatial directions, and anti-periodic boundary 
conditions, consistent with the imaginary time formalism used in (Q), in the temporal 
direction. 
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Figure 15: Plot of the fermion space-slice propagator, for 5.00 < k < 10.00 on a 16^ lattice, 
along with the fitted curves. Note how the sawtooth behaviour becomes more prominent 
as At increases, except on timeslices between 6 and 10, where noise dominates. The errors 
here are unbinned. 



Figure Uni shows the absolute values of propagators in the y-direction for k > 5.0. 
The propagators exhibit a more pronounced form of the sawtooth behaviour than 
Cfr [Cfx do not exhibit this behaviour at all). Unlike in the case of the pions in 
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[T3] . this is not due to the fermion becoming hght, as one can surmise from the slope 
of the curve. 



K 


TO 


X^d.o.f. 


fit window 


1.00 


0.97(1) 


2.540 


1-15 


3.00 


1.6(2) 


1.839 


2-14 


4.00 


2.1(1) 


1.233 


1-15 


5.00 


2.7(2) 


0.579 


1-15 


6.00 


2.6(2) 


1.097 


1-15 


7.00 


3.6(7) 


0.800 


1-15 


10.00 


5(2) 


0.793 


1-15 



Table 6: Fermion masses mfx in the x direction on a 16'^ lattice, with m = 0.03. 




Figure 16: Fermion screening masses ruf^, ''^fy^ ^.nd the geometric mean ^Jmf^mJ^^ 
versus k, on a 16^ lattice. Lines do not correspond to fits. 

As for the sawtoothed propagators in the r-direction, we performed fits only to 
odd ?/, as four-parameter fits proved unstable. The resulting screening masses are 
shown in tables IHl and [3 We can see that the fermions follow the same general trend 
as K increases as the pions in ^ those in the x-direction grow heavier, and those 
in the y-direction grow lighter; anomalies in the data (e.g. at k = 6.00) are likely to 
be due to noise. 

The geometrical mean y/ruf xfrif y increases from 1 to ~ 1.75(20) as we move into 
regions of large anisotropy, which suggests that some small dynamical effect may come 
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K 


m 


X fd.o.f. 


lit window 


nt 


i.UU 


1.1(2) 


U.oo 1 


A 1 O 

4-iz 


all timeslices 


3.00 


1.04(5) 


0.951 


3-13 




4.00 


0.9(1) 


0.907 


4-12 




5.00 


0.82(1) 


2.475 


1-15 


fit only 


6.00 


0.80(1) 


1.794 


1-15 


odd timeslices 


7.00 


0.80(1) 


0.689 


1-15 




10.00 


0.80(1) 


0.405 


1-15 



Table 7: Fermion masses 



in the y direction on a 16^ lattice, with 



m 



0.03. 



into play over and above that of the anisotropies themselves. This could correspond 
to a renormalisation of the parameter 6 in Lee and Herbut's model to a value other 
than unity, as a/ {6Xa){6X~^a) = 6a . 

5.2.3 Renormalised anisotropy: 




Figure 17: The renormalised together with a fitted curve, on a 16^ lattice. 



The renormalised fermion anisotropy HRf = rufx/mfy is displayed in Figure IT7I 
It is a measurement of the relevance of k > 1 relative to the particle in question. We 
find that the anisotropy parameter for fermions is irrelevant in the renormalisation 
group sense; that is 



R 



f 



1) 



< 1. 



(35) 
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Indeed, if we fit the above function to the data for 1 < k, < 7.00, we find Rf = .41(4), 
with a x^/d-o.f. of 1.67. This is striking, as the anisotropy is quite clearly relevant 
in the cases of pions ^3], and scalars as shown in ^4.31 The implication is that 
as K increases, the fermion - anti-fermion bound states become increasingly 1 + 1- 
dimensional, only able to propagate in the y — t plane (in the original condensed 
matter-inspired model // excitations associated with the other "flavour" , ie. node 
pair, would be confined to the x — t plane). The only excitations able to explore the 
whole 2 + 1-dimensional space are the elementary fermions. This point will be further 
discussed below. 

6 Discussion 

Here we summarise the main results of our study, and speculate as to the behaviour 
of QED3 as the anisotropy k is increased. 

We applied a finite volume scaling analysis to data from 16^, 20^ and 24'^ systems 
in an attempt to determine the order of the phase transition. There is no evidence 
for a diverging susceptibility as the volume increases, and remarkably the value 
marking the apparent transition appears to be very sensitive to system size; our 
fits assuming an isotropic model of finite volume corrections yield Kc = 7.66(5), 
which once the possibility of anisotropic corrections is admitted drifts out to Kc = 
12.3(6). Since the latter value lies outside our range of simulated parameters, it casts 
doubt on our original claim ^5] that a true chiral symmetry restoring transition is 
taking place. Rather, an interpretation of the transition in terms of a crossover from 
strong to weak coupling regimes seems admissible - very similar to the transition 
observed in simulations of isotropic QED3 j3 IIZ|- As in those studies, it appears 
to be a very difficult task to determine computationally whether chiral symmetry is 
actually broken in the weak coupling regime, reflecting the fact that QED3 may be a 
model with an abnormally large separation between the scale of dynamical symmetry 
breaking S and the natural mass scale g"^. It should be stressed, however, that the 
studies of the pion and scalar spectra in Sec. |3 are consistent with a chirally-restored 
vacuum at large n. 

What does seem clear is that any successful model of the finite volume scaling 
must take anisotropy into account - here our analysis assumed weak anisotropy, but 
models with differing critical exponents in different directions cannot be excluded. 
Unfortunately, the cure for these many uncertainties is to accumulate data from 
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many more values of L and m, which is beyond our current resources. 

However, it is intriguing to note that from fSllS], we can estimate that at T = 
we enter the dSC phase (and QED3 ceases to be a vahd effective field theory of the 
cuprates) somewhere in the region 6 ^ k ^ 8. Even taking the isotropic estimate 
Kc = 7.66(5) as the correct one, therefore, it is uncertain whether the intermediate 
pseudogap phase between SDW and dSC can actually exist. This raises a matter 
of some importance to future research: if k affects the behaviour of the system^, 
does it do so enough to make a difference in the condensed matter systems for which 
anisotropic QED3 is intended as an effective theory? 

Our studies of the propagation of // bound states in the scalar channel showed 
evidence for degeneracy between scalar and pseudoscalar as k increases, although 
the propagator data are markedly noisier in the scalar case. This is consistent with 
chiral symmetry restoration, but bearing in mind the cautious note of the preceding 
paragraphs, we should note that a very soft symmetry breaking cannot be excluded. 
Another important result is that the renormalised anisotropy kro- ~ urt^, implying 
that anisotropy is a relevant perturbation for both sets of particles. More graphically, 
this means that for large k ff bound states are effectively constrained to propagate 
in just the y-direction, and their dynamics are essentially 1+1 dimensional. 

The most significant result has emerged in the fermion sector, where we have found 
evidence that dynamical mass generation persists even once the apparent restoration 
of chiral symmetry has set in. Note that this result explains a rather surprising 
result reported in jT3]; namely, the average plaquette action f (O^,^) increases with 
K, implying that screening due to virtual // pairs in the quantum vacuum actually 
decreases with k, in contradiction to what would be expected if light fermion degrees 
of freedom were important in the high-K regime. The sawtooth structure that develops 
as K, increases may also be a sign of chiral symmetry restoration, although a study with 
m varying, beyond our current resources, would be needed to confirm this hypothesis. 

The fact that a non-zero dynamically generated fermion mass accompanies the 
chirally restored phase suggests that the symmetric phase is of an unusual kind. 
Witten |1T] has examined a similar situation in the Gross-Neveu model in 1 + 1 
dimensions; (xx) = and a dynamically generated mass may coexist if the following 
are the case: 

• The physical fermion is a branch-cut, not a pole, in momentum space, and lacks 
the same quantum numbers as the bare, massless fermion field ijj (the former 
■^That is, if k is relevant, or if k is irrelevant but Nfc is not universal. 
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has zero chirality - i.e. is chirally neutral - whereas the bare fermion has a non- 
zero chirahty). It follows that chiral symmetry tells us nothing about the value 
of the dynamically generated fermion mass - the system behaves in a chirally 
symmetric fashion in most respects apart from the existence of this mass. 

• There exists a massless (pseudo) scalar boson, which interacts strongly with the 
fermion and carries the chiral current. Interactions between it and the ip field 
(which is distinct from the observed physical fermion) are chirality changing. 
The interaction between ip and the boson modifies the chirally asymmetric por- 
tion of the fermion propagator, causing it to vanish. 

It is important to note that the scalar field in this example is not a Goldstone 
boson, which must be weakly-interacting. In the 1 + 1 dimensional Gross-Neveu 
model the formation of Goldstone bosons is prohibited by the Coleman- Mermin- 
Wagner theorem HH] , which states the impossibility of spontaneously breaking a 
continuous global symmetry in 1 + 1 dimensions. A similar phenomenon has been 
observed in simulations of the 2+ld Gross-Neveu model at non-zero T "H]. While 
perhaps it not clear how to define the effective dimensionality of an anisotropic theory, 
we take from this analogy the notion that infra-red fluctuations remain important in 
the chirally symmetric phase; in other words the interaction between fermion and 
scalar degrees of freedom is strong. 

In support of this hypothesis applying in the current situation, we point out 
that the mass ratios ruf : itIt^ : mo- vary from 1 : 0.2 : 0.4 at k = 1, consistent 
with broken chiral symmetry, to 1.8 : 0.9 : 0.9 at k = 10, consistent with restored 
chiral symmetry, but in which the scalar bound states are still tightly bound and 
light compared with the fermion mass scale. This should be contrasted with the 
"orthodox" chiral symmetry restored scenario m/ : m-,^ : ma ~ 0.5 : 1 : 1 observed in 
the 3(i Thirring model and portrayed in Fig. 17 of |21]. 

It should be noted that the situation in which dynamical mass generation without 
symmetry breaking is observed is sometimes referred to as pseudogap behaviour^. We 
must caution against confusing this with the pseudogap phase of the cuprate which we 
are modelling; while they share some behaviour in common (in both cases, we observe 
the phase disordering of an order parameter), they refer to different phenomena - the 
former referring to a phase of the putative effective theory, and the latter to that of 
the behaviour of the full description of the superconductor from which it is derived. 
^We thank Kurt Langfeld for bringing this to our attention. 
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Another important observation in the fermion sector is that K^f < k, imply- 
ing anisotropy is an irrelevant perturbation, i.e. fermions remain 2+ld particles 
as K increases, although scalar-mediated interactions among the fermions must be 
anisotropic. It will be a theoretical challenge to formulate an effective description 
incorporating these features. 

In many ways our study has raised more questions than it has answered; its main 
results have not been predicted by analytic treatments of the system performed so far. 
This may raise questions regarding the conception of the pseudogap in those models 
of HTc superconductivity - such as that of ^ - which require the presence of massless 
fermions in the chirally symmetric phase, since it appears that the expected link be- 
tween a non-vanishing chiral condensate and a dynamically generated fermion mass 
is broken. However, it is too early to make definitive statements; the fermion propa- 
gator should be measured in a number of gauges, so that we can be certain as to how 
much (if any) of the observed behaviour is an artefact of Landau gauge. Ultimately, 
more data on how the dynamically generated fermion mass behaves as the chiral, 
thermodynamic and continuum limits are approached will be needed. Anisotropic 
QED3 appears to be every bit as computationally demanding and as fascinating as 
its isotropic counterpart. 
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